Negative Entropy Production in Oscillatory Processes 
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Linear irreversible thermodynamics asserts that the instantaneous local spontaneous entropy pro- 
duction is always nonnegative. However for a viscoelastic fluid this is not always the case. Given 
the fundamental status of the Second Law, this presents a problem. We provide a new derivation 
of the Second Law, from first principles, which is valid for the appropriately time averaged en- 
tropy production allowing the instantaneous entropy production to be negative for short intervals 
O ■ of time. We show that time averages (rather than instantaneous values) of the entropy production 

' are nonnegative. We illustrate this using molecular dynamics simulations of oscillatory shear. 

Linear irreversible thermodynamics asserts that in local equilibrium p[, the spontaneous entropy production per 
unit time, per unit volume, the so-called entropy source strength cr(r, t), cannot be negative @, 0]. Further it states 
that the entropy source strength is a sum of products of irreversible thermodynamic fluxes Ji and forces Xi , 
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a(r,t)=^J i (r,t)X i (r,t)>0 W 

where the source strength is calculated at a position r and a time t. For steady state processes close to equilibrium 
Eq. [T] is clearly correct. However, for processes that involve time dependent or oscillatory thermodynamic forces, in 
viscoelastic materials close to equilibrium, Eq. [T]is incorrect. These problems are usually resolved by separating the 
fluxes into "storage and loss" components Q. Such a separation is not derived from first principles and is process 
specific. A second problem with Eq. [TJis that it is restricted to the near equilibrium, linear response, regime. Further 
from equilibrium the definition of entropy and therefore temperature remains an unsolved problem [5(. Until very 
recently there has been no known generalization of Eq. [T]to the far from equilibrium regime. 

Recent advances, which are of broad interest @, 0|, have cast new light on these issues. In particular we will use 
the Jarzynski Equality (JE) 0, [§] to derive a variation on Eq. Q] from first principles. In order to proceed we first 
define the conjugate flux J, for a system which is driven away from equilibrium by a dissipative field F e , in terms of 
the rate that work is done on the system [Io| . 

> ■ 

^ JV-F e (t) = -k B Tn(t)=k B TA(T(t))-j t H (T(t)). (2) 

O ! 

^vq , Here V is the system volume, k B is Boltzmann's constant, T is the temperature of the synthetic thermostat or a 
j— ( ■ large heat reservoir [ll[, fi(t) is the dissipation function as defined for the Evans-Searles Fluctuation Theorem fl2|, 
r(t) is the phase space vector Hq is the internal energy, k B TA gives the rate that heat is lost from the system to 
the thermostat or reservoir [ll[ and A is the phase space compression factor [hJ • Eq. [2] is defined for arbitrary field 
' strengths and in the linear response regime it gives the entropy production, 
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(JV ■ F e (t))/T = k B (fi(t)> = / dr a (r, t) + O (F e 4 ) , (3) 
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where (...) denotes an ensemble average. 

The JE has been shown d @] and proved El for the time reversible thermostatted dynamics we employ below. 
It has also been experimentally verified [LJ, LL5j. The JE gives the change in Helmholtz free energy AA for a 
system which has undergone a nonequilibrium process starting from an initial (t — 0) equilibrium distribution of 
phases /i(r,0) oc exp[— j3H(T, A(0))] and initial Hamiltonian H(F, A(0)) to a final Hamiltonian H(T, A(r)) at t = r 
and then relaxed to a new equilibrium. The dynamical proof [Lj| requires that the dynamical pathway F(t) be 
such that the distribution of phases at t = r can subsequently relax to a final canonical distribution at t = oo, 
/oo(r, oo) oc exp[—/3H(T,X(T))]. The parametric transformation of the Hamiltonian is complete at a finite time 
t = t, by which time the system is not expected to have fully relaxed to the new equilibrium. JE states, 
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(exp (-0W T )) = cxp (-0AA) . (4) 

where W T = J Q T dsW(s) is the Jarzynski work function and W(t) = [^H (T(t), A(t))] - fc B TA(r(t)), Although 
H (T(t), A(t)) ^ H (T(oo\ X(t)) it is clear that, W T = W^. 

The proof of the JE [l3| requires that the two systems be connected by a path 1 — > 2 and its inverse path 2 — > 1 . 
The proof does not put restrictions on the time dependence of the path. The parametric change in the Hamiltonian 
from A(0) — > A(r) may in addition contain work due to the system being driven by a dissipative external field [T3|. If 
the work is solely due to a dissipative external field (A = V t) then the rate of work W will be the same as that 
given by Eq. [5] (ie W = k B TQ), 

The JE allows a first principles proof of the Second Law Inequality (SLI) from the equations of motion as first 
shown by Jarzynski 0|. By combining Eq. [4] with the mathematical identity exp(cc) 1 + x we have, 

e -l3AA = ( e -P(Wr-{W T ))j e -P{W T ) > g-,8<W,-> _ ^ 

Noting that is a monotonically increasing function we derive the Clausius Inequality 

(W T ) > AA, V t > 0. (6) 

If the system of interest is a fluid and if the transformation involves, say, a shearing deformation or perhaps the 
translation of one particle through a fixed distance, then clearly AA — 0, we may treat the field as external A = 0, 
and thus, by Eq. W T = k B TVL T T = - ds J (s) V ■ F e (s) where Cl T = \ Q ds fi(s). The SLI then follows, 



(ti T ) = -i- / ds (J(s)V -F e (s)} >0, V r>0, V F e (s) (7) 

which forms a generalization of Eq. [TJ that is valid at arbitrary field strengths and derived from the equations of 
motion. A significant difference between Eq. [7]& Eq. [T]is that Eq. Q] applies to instantaneous values whereas Eq. [7| 
applies to time averages of the entropy production starting at t = from an initial canonical distribution. We note 
that Eq. [3 ha s previously been derived under a more restrictive set of conditions from the Evans Searles Fluctuation 
Theorem |lGl |. 

We decided to test this prediction using nonequilibrium molecular dynamics simulations of shear flow in a fluid. 
We consider the case of sinusoidal shear applied to a viscoelastic fluid. We employ the Lees-Edwards (sliding brick) 
periodic boundary conditions along with the so-called SLLOD equations of motion for planar Couette flow [id ]. 

Qi = — + H(t)q yt (8) 
m 

Pi = Fj - i-y {t)p yi - apt 

where j = du x /dy is the strain rate, Pi is the peculiar momentum taken relative to the streaming velocity u x (y) = 
ij(t)y and a(t) is a Gaussian thermostat multiplier which holds the kinetic temperature fixed [[T3] - For Couette flow 
with a constant strain rate the SLLOD equations of motion are known to give an exact description of steady adiabatic 
planar Couette flow for arbitrary values of the strain rate [lfj. For time dependent shear flows it is known that 
the SLLOD equations give an exact description of such flows in the linear response regime for both adiabatic and 
thermostatted flows [XQ] . For high shear, oscillatory flows it is not known whether SLLOD is exact but it is widely 
assumed this is so. 

To apply the SLI, Eq. \7\ to the equations of motion Eq. we note that the equilibrium Hamiltonian or internal 
energy is Hq = Y^Pi/^ m + $ where $ is the total interparticle pair potential. The dissipation function, which in the 
linear response regime gives the entropy production Eq. may then be obtained from, 



N N 

k B Til (t) = -jVP X y - -7(E V -^T ~ E F *iJ9vij) . (9) 

i—l i<j 

where N is the number of particles, P xy is the xy element of the pressure tensor, F xi j is the x component of the 
pairwise additive force on particle i due to particle j, q y ij is the y component of the vector connecting their centers 
and p X i is the x component of the momentum of particle i. So in this case for JV ■ F e we have J = P xy and F e = 7. 
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Figure 1: The response of the system to the oscillatory shear, in the linear response regime, which began at time t = following 
Eq. The solid line is the entropy production, = — jP xy (t)V/e. The dashed line is the instantaneous rate of heat 

absorbed by the thermostat (dQ/dt) /e. The dotted line gives the strain rate j(t) multiplied by a factor of 10. The data for 
the heat exchange is considerably more noisy than the other data. Clearly (fi(t)) is at times negative. 
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Figure 2: The dashed line is the integral, (fit) t, which is nonmonotonically increasing but always nonnegative. The solid line 
(black) is the time average (fit) which is also always nonnegative. 

Simulations of oscillatory shear were carried out on a fluid in three dimensions using the pair potential 4>ij — 
e [( a / r ij) 12 ~ ( <J / r ij) 6 + 0.25] V nj < 2 1 / 6 , with the volume and temperature held constant. The equations of motion 
were solved using a fourth order Runge-Kutta algorithm with a tim e step At = 0.001. The number density is 
p = Na 3 /V = 0.95, N — 108, ksT = e and the time unit is y^ma 2 Je throughout. For times t < the system was 
in equilibrium. The initial equilibrium configurations were obtained by sampling an equilibrium trajectory, at time 
intervals of 5, which was generated by solving Eq. [8] with j = 0. Starting from these initial configurations a total of 
5 x 10 5 nonequilibrium oscillatory Couette flow trajectories were computed: j(t) — 7osin(u;i), t > with id = 4tt 
and 70 = 0.2. The duration of the nonequilibrium trajectories was 2. 

In Fig. 1 the dissipation function, (f2(£)) m units of fee (which was calculated by ensemble averaging all of the 
nonequilibrium trajectories) is plotted as a function of time along with the ensemble averaged rate of heat absorbed 
by the thermostat, Q(t) = —3NkBT(a(t)) = ksT (A) in units of e. Also shown is the strain rate, 10 x j(t). The 
initial transients in the response decay very rapidly. It is clear that the response of the fluid is viscoelastic: there is 
a phase lag between P X y(t) and j(t) due to the relatively high frequency u>. This shows that the dissipation function 
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(or at weak fields, equivalently, the entropy production Eq. [3]) is negative within certain intervals of time even though 
the system's response is linear. This effect is not due to the amplitude of the shear rate 70 being large. This is clearly 
at odds with the traditional view from_ irreversible thermodynamics. In contrast Eq. [7] is satisfied at all times as may 
be seen in Fig. 2 where the integral (fit) t and the time average (fit) are plotted as a function of time t. 

In summary we have shown the assertion of linear irreversible thermodynamics that the instantaneous entropy 
production is always nonnegative is incorrect for the case of time dependent viscoelastic fluids even if they are in 
the linear response regime close to equilibrium. The Second Law Inequality Eq. [7] derived from the Jarzynski 
Equality states that the time integral (starting from t = 0) of the ensemble averaged dissipation function cannot be 
negative for arbitrary integration times and arbitrary field strengths (of course in the weak field limit the dissipation 
function is equal to the entropy production Eq. [3]) . This inequality requires that the Helmholtz free energy of the 
corresponding equilibrium system does not change. For planar shear this is a necessary condition for the fluid state, 
which by definition cannot support a constant stress. For a solid, undergoing oscillatory nonplastic deformation, the 
equilibrium free energy would depend on the deformation and the Clausius Inequality Eq. [6] will need to be used 
rather than the Second Law Inequality Eq. 

Lastly we note that the Second Law Inequality is a macroscopic consequence of the Jarzynski Equality and of the 
Evans Searles Fluctuation Theorem [16]. All previously derived consequences of the JE and the Fluctuation Theorem 
were microscopic in nature. The Second Law Inequality in the form Eq. ]7\ has important consequences in applications 
such as atmospheric physics where the principle of maximum entropy in nonequilibrium states has been employed 
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